Two-phonon spin-lattice relaxation of rigid atomic clusters 
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Spin-phonon relaxation due to two-phonon processes of a spin cluster embedded in an elastic 
medium has been studied. For the case of uniaxial anisotropy, relaxation rates due to Raman 
processes and processes involving the emission of two phonons have been obtained. For a biaxial 
spin Hamiltonian, the rates of transitions between tunnel-split levels have been computed. By 
comparison with the rates of corresponding direct processes, we have established temperature ranges 
where the Raman mechanism dominates over the one-phonon decay mechanism. 
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I. INTRODUCTION 

Relaxation of spins in a paramagnetic solid is a prob- 
lem of fundamental interest. It is also related to impor- 
tant applications such as spin resonance and the use of 
spin as a qubit. In a paramagnet, the decay of a spin state 
is due to the interactions with phonons, nuclear spins, 
dipolar fields, etc. In principle, one could suppress the 
interaction with nuclear spins and dipolar fields but the 
interaction with phonons always remains. Often it domi- 
nates spin relaxation. Thus, spin-lattice interactions pro- 
vide the most fundamental upper bound for the lifetime 
of spin states in paramagnets. The spin relaxation due to 
interaction with the lattice can occur by various mecha- 
nisms. The most studied and often the dominant ones are 
direct processes, in which a single quantum is exchanged 
between the spin system and the lattice. It was pointed 
out by Waller in studying the modulation of the spin-spin 
interaction by the lattice waves that, unlike in electro- 
magnetic phenomena, the inelastic scattering of a phonon 
combined with a transition in the spin system could be 
very important jl]. This is a two-phonon process con- 
sisting of the absorption of one phonon and the emission 
of another phonon with different frequency. The mech- 
anism is analogous to the Raman effect in optical spec- 
troscopy and it is often referred to as a Raman process. 
Spin-lattice relaxation mechanisms based on dipolar in- 
teractions (Waller) were insufficient, however, to account 
for the transition rates measured in experiment. Hcitlcr 
and Teller considered Raman processes based on a 
more potent mechanism based upon modulation of the 
crystal electric field under the action of the lattice vibra- 
tions. Their theory was further developed by Kronig 3] 
and Van Vleck .4]. They obtained a spin-phonon cou- 
pling based on the spin-orbit interaction that permitted 
calculation of relaxation rates that were of the same order 
of magnitude as the experimental ones. Later on, Orbach 
P| managed to simplify the treatment of the problem by 
expanding the crystal electric potential in powers of the 
fluctuating strain caused by the lattice vibrations. 

Recently, this problem has received new attention in 
connection with spin relaxation of molecular clusters. 
The spin of many such clusters is formed inside a rel- 



atively rigid magnetic core that can rotate in the pres- 
ence of the deformation field but is more resistant to 
distortions of the core itself. It has been noticed that the 
spin relaxation of such a cluster can be obtained within 
a model that is parameter free, that is, it gives the relax- 
ation rates in terms of the known crystal field Hamilto- 
nian of the magnetic core Ha |g] ■ Even for non-rigid clus- 
ters, calculation of the effect of rotations is meaningful. 
It has been theoretically established that spin-phonon re- 
laxation rates due to both, one-phonon and multi-phonon 
processes, are inversely proportional to some high pow- 
ers of the sound velocity [||. Since longitudinal phonons 
have a larger sound velocity than the transverse phonons, 
processes involving longitudinal phonons can be safely 
neglected. The effects of the transverse phonons can be 
split into shear deformations of the lattice cell and local 
rotations of the lattice that preserve the symmetry of the 
crystal field. To describe deformations of the first kind 
one needs to employ terms in the Hamiltonian contain- 
ing phenomenological coupling constants, whereas the lo- 
cal rotations can be described by a parameter-free spin- 
phonon Hamiltonian that is defined solely by the form of 
Ha- In general, processes due to the shear distortion of 
the lattice and those due to the local rotation of the lat- 
tice should result in comparable relaxation rates. Even in 
this case, the latter are of a fundamental importance be- 
cause they provide a parameter-free lower bound on the 
decoherence of any spin-based qubit. For rigid spin clus- 
ters, interaction of the spin with rotations of the crystal 
field is the only source of spin-lattice relaxation. 

The angle of rotation of the crystal field axes 
(eW, e^ 2 \ e^ 3 )) due to the action of the phonon field u(r) 
is given by 



S(f>(r) = — V x u(r). 



(1) 



Because Ha is a scalar, the rotation of (e^, e^ 2 ), e^ 3 )) is 
equivalent to the rotation of the vector S in the opposite 
direction. As it is known, the rotation of the operator S 
can be performed by the (2S + 1) x (2S + 1) matrix in 
the spin space, 



RSR- 



R = 



-iS-S<p 



(2) 



2 



Then, the total Hamiltonian in the presence of phonons 
can be written in the form |8l 



H = RH A R~ 1 +H Z +H 



ph 



(3) 



where Ha is the crystal-field Hamiltonian in the absence 
of phonons, Hz is the Zeeman Hamitonian and H p h is 
the Hamiltonian of harmonic phonons. In these formu- 
las 8<p must be treated as an operator. The canonical 
quantization of phonons and Eq. Q yield 



8^ 



E 

kA 



\ik x eic\l e lk ' ; 



(a k A + a f _ kA ) (4) 



where p is the mass density, V is the volume of the 
crystal, are unit polarization vectors, A = t\,ti,l 
denotes polarization, and lu^x = v\k is the phonon 
frequency. 

For one-phonon processes, relaxation rates of spin- 
phonon transitions have been recently computed with 
the help of the above formalism in Ref. 0. Such pro- 
cesses dominate spin-phonon relaxation at zero tempera- 
ture when no thermal phonons are present in the system. 
At finite temperature, however, two-phonon processes 
may take over. In this paper we study two kinds of two- 
phonon processes. The first kind consists of an inelastic 
scattering of phonons by the spin-system, or spin-phonon 
Raman processes. It corresponds to the annihilation of 
an incoming phonon of frequency Wk and the creation 
of an outgoing phonon of frequency w q , with k and q 
being the corresponding wave vectors. The second kind 
involves emission of two phonons. Note that the conser- 
vation of the energy in spin-phonon interactions requires 



hv„ ± ftwk = AE > , 



(5) 



where AE is the energy difference between the spin-states 
and wi(,Wq are the frequencies of the phonons involved 
in the transition. The plus sign applies to processes in- 
volving the emission of two phonons and the minus sign 
applies to Raman processes. When AE <C fiuo (with 
u>d being the Debye frequency), Eq. JSJ causes the phase 
space of phonons to be much greater for Raman processes 
than for processes involving the emission of two phonons. 
Consequently, the spin-phonon Raman scattering usually 
dominates over the processes involving the emission of 
two phonons. However, the same condition JSJ implies 
that the energy of the phonon emitted in the Raman scat- 
tering process, hu>q, must satisfy hujq > AE, whereas 
in the process involving the creation of two phonons 
their energy must be smaller than AE. Therefore, if 
fcgT <C AE, the total number of phonons available to 
carry out the spin transition may be much greater in 
the two-phonon emission case than in the Raman case, 
so that the former case can become dominant. In both 
cases the matrix element of the transition is a sum of 
two terms. The first term, M^- 2 ' , comes from the first 



order of the perturbation theory on the spin-phonon cou- 
pling containing a product of two phonon displacement 
fields. The second term, M^ 1+1 \ comes from the sec- 
ond order of the perturbation theory on the spin-phonon 
coupling that is linear in the phonon displacement field. 
In some cases these two terms interfere, so that the re- 
sulting transition rate, based upon \M^ + A/( 1+1 '| 2 , is 
different from the one obtained by adding up the rates, 
|M( 2 )| 2 and |M^ 1+1 )| 2 , that each term would produce by 
itself, as was incorrectly done in the past 



II. MATRIX ELEMENTS OF TWO-PHONON 
PROCESSES 

The treatment of two-phonon processes requires con- 
sideration of terms up to second order in phonon ampli- 
tudes in the Hamiltonian: 



RHaR 1 — Ha + Hs-ph 
with 



H s -ph — H s 



(i) 



- P h 



H 



n 



(1) 

s— ph 



(2) 

s— ph 



2! 



Ha, S a 



:S, 



The total Hamiltonian can then be written as 
H = Ho + H s -ph , 



(2) 

s — ph ' 

(6) 



(7) 



(8) 



where Ho is the Hamiltonian of non-interacting spin and 
phonons, 



and 



Hs — Ha + Hi 



(9) 



(10) 



is the spin Hamiltonian. 

We will study spin-phonon transitions between the eigen- 
states of Ho, which are direct products of the spin and 
phonon states, 



k±) ® \4>±) 



(ii) 



Here \tp±) are the eigenstatcs of Hs with eigenvalues E± 

(E+ > E-) and \<t>±) are the eigenstates of 7i p h with 
energies -E p h,±- Spin-phonon transitions conserve energy 



E.+E, 



P h, 



E_+E. 



ph, 



(12) 



To obtain the relaxation rate of the transition |^+) — > 
|^-) one needs to evaluate the matrix element of the 
process. This matrix element is the sum of the matrix 
element with H^} ph and that with H^l ph in the second 
order Hal: 



M = M (2) + M (1+1) 



(13) 



3 



where 



and 



M (2) = ( * 



n. 



(2) 
s — ph 



(14) 



(i) 

s— ph 



H 



(i) 

s— ph 



(15) 



Here £ labels intermediate spin-phonon states. 

III. RAMAN PROCESSES 

For the Raman processes of interest, a phonon with the 
wave vector k is absorbed and a phonon with the wave 
vector q is emitted. We will use the following designa- 
tions for the phonon states 

\<i>+) = l^ki^q) , \<f>-) = \rik - l,n q + 1) (16) 

In this case, the conservation of the energy reads: 

AE = E + - E- = hujq - hoj k . (17) 

The matrix element of the Raman process can be writ- 
ten as Mr = M R 2) + M R 1+1) . According to equations Q 
and UU) 



Ha, S a 

x (n k - l,n q + 1 1^(5^1 n k ,n q ) . (18) 
It is convenient to express the phonon matrix element as 



(nk - l,n q + 1 \54> a 
where 



"k 



» «q) = M ph-R 



M, 



/3a 

ph-R' 

(19) 



Kh-R = N + 1 "q> ( n * - 1 1^/9 1 »k) • (20) 
With the help of Eq. one obtains 



M, 



a/3 _ 
ph-R 



H 2 [kxe kAk ] Q [qxe qAq ] / 

^ftW k A k ftW q A q 



\J(riq + 1) n k . 



(21) 

On the other hand, using the definition l|2()(l . one obtains 
from Eq. (j*H)|l 



71/ 



(i+i) _ 



E 
E 



Ha, S a 



Ha, S/s 



E + + fit<j k — -E^ 



M, 



a/3 

ph-R 



Ha, s a 



Ha, S/3 



E + — Ef — ftlOn 



M. 



/3a 

ph-R ' 



(22) 

where £ now labels intermediate states of the spin only. 
The intermediate phonon states are |n k — l,^ q ) in the 
first term and |n k ,n q + 1) in the second term. 



A. Transition between eigenstates of S z 



Consider for example the spin Hamiltonian 
Us = H A + H z = -DS 2 Z - gfi B H z S z 



(23) 



that commutes with S z . The exact energy states of 
this Hamiltonian are the eigenstates of the S z operator, 
S z \m) = to | to), with energies given by 

E m = -Dm 2 - gfi B H z m . (24) 

Let us study the general case of the spin transition to + — > 
TYi— . From equations i|18|) . i|19[) an d (|23|l one obtains 



(25) 



a/3 



(rn + + m_ — 2m|)(TO_|S , Q ,|TO5)(TO5|S , ( 3|m+) 



/ 7 / 



with M 



i 

a/3 



Mpj^_ R ) . The summation 



i f M Q/3 

X ph-R ~~ 2 ^"ph-R 

on to.£ runs over all the spin states. On the other hand, 
Eq. results in 



M 



D 2 {m 2 _ — rn^){m\ — to|) 



(i+i) = 

'' ^ £>(to| — to+) + g/i B H(m^ — m + ) + fiw k 

x (m_ 1 5a | m ? ) (to c I S f3 1 m+) M^f_ R 



E 



D 2 {m 2 __ — to?)(to 2 k — TO?) 



D(m^ — m + ) + g(i B H(rri£ 



X (m_ | S a |m £ ) (m e 1 5/3 1 m+)M^_ R , (26) 
where summation over a and /3 is assumed. 

1. Adjacent spin levels, m — > ?n ± 1 

We will first treat the Raman processes involving tran- 
sitions between adjacent levels of the spin-Hamiltonian 
(|23|l (see Fig. Therefore, the spin-states in this case 
will be 

\m + ) = \m) , |m_) = \m ± 1) , (27) 

where the plus sign applies to positive to and the minus 
sign applies to negative to. A straightforward calculation 
of the matrix elements in equations (|25|) and l|26|l leads 
to 



M 



±z 



M R 1+1) = , (28) 



where 



(29) 



and l m , m ±i = \/{S Tm)(S ± to + 1). 

(2) 

Then Mr = M R and the transition rate is given by 



4 



w -10 




w -10 



-5-3-11 3 5 




FIG. 1: Transition between adjacent spin-energy levels of 
Hamiltonian I23H mediated by Raman processes. 



FIG. 2: Transition between the spin-energy levels of Hamil- 
tonian 1231 1 m — > m + 2 due to Raman scattering. 



pm-miX = ^ ^|M^| 2 «5 (fiw q - ft Wk + E m±1 - Em) . 

kA k 
qA q 

(30) 

Note that in the sums over the polarizations A k and A q , 
only the two transverse modes are considered. To com- 
plete the calculation, we make use of 

[k x eict] = ±fce k f (31) 

E(k • a)(k • b) . . 
(e kt -a)(e kf • b) = a b-^ g '- (32) 

t=ti,i 2 

and the replacement of ^ k by V J d 3 k/(2ir) 3 to obtain 

r m^m±i 1 4,m±i / m [D(-2m T I)] 2 ^ 7 
Tr = h—^ £f {kBT) ' (33) 



R 

where 



£t = pvt) 



3 5\l/4 



(34) 



is a characteristic energy in the problem. In these expres- 
sions p is the mass density and Vt is velocity of transverse 
sound. Iri is given by 



Im = 



1 



1152 Jo 



e D /T 



dxx 3 (x + ei) 3 



,'>' + <" 1 



where 



{e x - l)(e x+£ i - 1) ' 
(35) 



E m -E m±1 D{±2m+l)±gii B H 
£ i = ; — ~pp, = j—^ ■ ( 36 ) 



k R T 



We remind the reader that in the above formulas the 
choice of upper and lower sign corresponds the choice of 
± in Eq. (|27|l . We use to + 1 for positive to and we use 
to — 1 for negative to, so that E m — E m ±i > 0. 



2. Non-adjacent spin levels, m — > m ± 2 

It is clear from equations (|25|l and 1)26(1 that the only 
allowed transitions between non-adjacent spin levels are 



those with 



\m+) = \m) 



\m±2) 



(37) 



In this case, represented in Fig. [21 equations and 
(|2FJ|l give 



M 



(i+i) 



D' 2 



(±2to + 3)(±2to -I- 1) i m ,m±l^m±l,m±2 



M. 



±± 



ph-R. 



M 



ph-R 



^m±l — -Em + fi-^k ^m±l — E ri 



(38) 



where 










M ph-R 


1 

2 




H^ p \ X -r) 


JW ph-R 


^pV-R 


1 

2 




H ^h-R.) 


- M vv 



(39) 



The transition rate can be obtained by computing 



-\m — >m±2 
R 



2tt 

^ -t-|M_r| 2 <5 (7iu; q - Tiw^ + E m ± 2 - E m ) , 



kA k 



(40) 

with Mr given by equations (|13[) and (|38fl . Again, wc 
use equations lj3*T)l and ijS*^ and the replacement of ^ k 
by Vd 3 fc/(27r) 3 (with V being the volume of the crystal) 
to obtain the final result 



— >m±2 



where 



1 ^m,m±l^m±l,m±2 D 2 
H 7T 3 £? 



(A: B T) 7 , (41) 



1 /" 9d/T , (z + e 2 ) 3 a; 3 e( a;+£2 ) f 
fi2 ~ 288 7 dx [ e --l][e(^)-l]l 



1 



(±2to + 1)(±2to + 3) 



ei + x ei — £2 — a; 



(42) 



5 




where A is the tunnel splitting of the levels to and to' 
that can be calculated from the exact spin Hamiltonian 
Hs ^2 or determined experimentally. Diagonalizing this 
2x2 matrix yields the eigenvalues 



E± = ~ (E m + E m , ± ^Jw 2 + A 2 



(48) 



The corresponding eigenvectors can be represented in the 
form 



1 



|V> ± > = ^(C±|<U±CM4w>) 



(49) 



where 



FIG. 3: Transition between tunnel-split states due to Raman 
scattering. 



c± = Ji± 



w 



y/W 2 + A 2 



(50) 



with 



E m -E m±2 4:{±m+l)±2g(i B H 
e 2 — ; — — — ; — ^ • ) 



k B T 



k B T 



Far from the resonance, \W\ ^> A, the eigenstates and 
energy eigenvalues reduce to those of \<fi m ) and \<p m /) 
states. 



B. Transitions between tunnel-split states 



2. Matrix elements 



Consider now a biaxial spin Hamiltonian with strong 
uniaxial anisotropy 



Hs — Ha - SM^H ' S , Ha 



-DS Z + E[S X - S y ] , 
(44) 

where D >• E > and DS 3> gfi B H±, with H± being 
the transverse magnetic field. Consequently, Ha nearly 
commutes with S z . The energy levels of this Hamiltonian 
are approximately given by Eq. I|24|) . The two levels m 
and to' are in resonance for the values of the magnetic 



field H 
by 



(res) 



(m + m')D/ (g^ B ). The level bias is given 



W 



E m - E' m = 9t i B (H z - H { ;Xl) (m' 



(45) 



1. The two state model 

Due to the terms in Hs that do not commute with S z , 
the true eigenstates of Hs far from a resonance are given 
by expansions over the complete |to) basis: 



\4>n 



i" = -S 



(46) 



where 



1 and all the other coefficients are small. 



Hybridization of the states |0 m ) and \4> m ,) when they 
are close to resonance can be taken into account in the 
framework of the two-state model 



. ) = E mi , m^m, to' 



(47) 



Here we consider Raman processes involving spin tran- 
sitions between tunnel-split states (see Fig- EJ . That is, 
the spin eigenstates in Eq. Hll|) are given by Eq. H49fl . 
In order to compute the matrix element of the Raman 

(2) 

process, we can rewrite M R by adding and subtracting 
Hz in the spin matrix element of Eq. I|18|) 



Ha, S a 



1>- 



Hz, S a 



Hs, S a 



(51) 



Taking into account that \ip±) are the eigenstates of Hs 
and inserting the identity 1 = j^) (^| we can express 
the first term of the right hand side as 



i>_ \[[Hs,S a ] ,Sp 



(E- +E + - 2E e ) \S a \1>e) \S P \ 1> + ) . (52) 



Then, 



71/ 



(2) 



J2(E-+E+-2E ( ){^_\S a \^z) 



(v € i^iv+) -U 



Hz, S a 



M, 



a/3 

ph-R • 

(53) 



Following the same procedure, we can rewrite 
from Eq. 11' I'D as an expansion on powers of H: 



(i+i) 



M 



(i+i) 



= M ( B 1+1) (H°) + M R i+1) (H 1 ) + 0(H 2 ) (54) 



r(i+i) 



6 



with 



m r 1+1 \h°) = 

-J2( E -- - E +) it- \s a \ ^) 1^1 

£ 



a/3 

ph-R 



M ph-R 



£+ + — E% E + — E/: — htu q 



(55) 



and 



E 

£ 



Hs, Sp 



E 

£ 

E 

£ 

E 

£ 



5 G 



E + + fitJk — 



Hz,S a 



-E+ — E^ — hujq 



H s , S a 



Hz,Sp 



E + — E^ — tkjjq 



M, 



a/3 

ph-R 



M 



a/3 



ph-R 



M, 



/3a 

ph-R 



M. 



/3a 

ph-R' 

(56) 



H = Wi X) ® W[ 2) , (-5, -5 + 2, . . . , 5) G ?4 X) and 

(-5+1,-5 + 3, 5-1) e Tif 3 . Thus Eq. JHEJ 
yields a zero result. 

Let us consider now the terms with £ ^ ±. In this 
case, the difference between the energies of the doublet 
is much smaller than the energy distance to the other 
states, \E+ — £L| <C \E± — E%\, so that one can re- 
place E± with E and hujq with huk in the matrix el- 
ements. We consider the case of low temperature, when 
hojk <C \E — E$\ for thermal phonons. Then Mr can be 
simplihed to 



M R ~ 2 (fr Wk ) 2 M^_ a Ac fi 



(61) 



with 



^, r <*-W*«X*.IW^> , ((B , 

£ C 

where prime means that £ = ± have been excluded. Note 
that the quadratic dependence on w q in Eq. (|61() results 

(2) 

from cancellations between terms from M R and terms 



from M R . Consequently, the relaxation rate will have 
a different temperature dependence from the result that 
one would obtain if one added the rates stemming from 
M R ^ and M R +1 ^ independently. The rate is given by 



3. Transition rate for H — 
At H = one obtains 

M * = - E l^£> (^£ i^i <M 

£ 



JW ph-R X 



E-+E+ (EL - E 6 ) (Eg - E+) 
o #£ + 



(/3a) 
ph-R 



E_+E, 



-E 



E + + htuk — E^ 

(E_ - E e ) (E e - £ + ) 



E + — E% — /iWq 



(57) 



It is convenient to consider the terms with £ = ± and 
£ 7^ ± separately. The contribution from £ = ± is 

M R = -(£+-£_) ((V-|5a|V>->(V'-| 1 W+) 



■<V>_|£y*A + )<V' + IW + ) M p t_ R 



(58) 



Using the time-reversal symmetry, we obtain 

(V ± |S|^ ± H-(V+|Sh/+>*. (59) 



V~ -EE y |m « |25 ^ ~ 

k,A k q,A q 



1 V^M' + A* ) (fcsT)n 



18fi(27r) d ^ 



Z" 8 



10 



(63) 



where equations (|31|l and 1321) have been used and the 
continuum limit has been taken. The constant I w is 



/in = 



dx ■ 



x w e x 



'10 — i 2 

/o (e* - 1) 



r(ii)c(io). (64) 



For transitions between the lowest doublet m = — 5, m! = 
5 specified in Eq. (|49|) one can evaluate A a p by consider- 
ing also the first excited doublet m = — 5+ 1, m' = 5—1 



|V 1± ) = -j= (Ci± |0 m+1 ) ± C lT l^^i)) (65) 

with Ci± given by Eq. (J5DJ with A — * Ai and T4 7 — > W\. 
Hence the only non zero matrix elements are 



For the biaxial model with E > the states \ip±) are 
real. Then, 



E m — E m +i 



-A 



(66) 



(^±|S,|V±)=<V'±IW ± )=0. 



(60) 



On the other hand, ( , 0_|_|5j,|V'±) = be- 
cause of the factorization of the Hilbcrt space: 



where 



7 



Evaluation of these matrix elements yields 

S 2 A 2 

E A aP {A* a0 + A; a ) = {E _ s _ E _ g+ifw2 + A2 • 

(68) 

For the case under consideration W = and the transi- 
tion rate is, then 



r 



_ 80tt 7 S 2 {k B T) n 



no 



(69) 



where 



(76) 



Again, it is convenient to replace K^ 1 ^ in Eq. (|75() by 
the symmetrized version 



(77) 



The addition to the Raman matrix element due to H ^ 
is, then 



4- Transition rate in the presence of a magnetic field 

Here we are going to evaluate the contribution of 
the magnetic field H to the rate of the transition be- 
tween the tunnel-split ground-state levels. We consider 
magnetic fields with very small longitudinal component, 
gfi B H z ~ A. For longitudinal fields gfi B H z ^> A the Ra- 
man processes die out. We are also restricted to not very 
large transverse magnetic fields, gn B H± <C DS. In this 
case, it is sufficient to compute the lowest-order contri- 
bution of the magnetic field to the transition rate. This 
contribution can be essential because of the cancellation 
that occurs in the matrix element for the H = case, 
Eq. (|61(l . This cancellation leads to the oc T 11 de- 

pendence. As we shall see, there is no such cancellation 
in the field-dependent term, so that the result shows a 
T 7 -dependence. 

According to Eq. (|53[) . the linear order contribution of 

(2) 

H to M R is given by 



SM 



(2) 



1>- 



Hz, S a 



(70) 



The double commutator equals 



Hz, S a 



g\i B (HpS a - HyS-fSap) (71) 



so that 
8Mg> 



where 



5Ms (HpS az - H z 6 a p)M"£_ R (ip_ \S Z \ ip + ) 



VW 2 + A 2 2 



(72) 



(73) 



The phonon matrix element M^_ R is symmetric in a/3. 

(2) 

Thus, it is possible and convenient to replace K a p in Eq. 
1|72JI by the symmetrized tensor 



k (2) _ 1 / (2) 



K 



(2) 



(74) 

On the other hand, the linear order contribution of H 
to Af£ +1) is given by Eq. (JHSJ. By using Eq. (|71|l and 

assuming hiu^ <C £m+i — E m we obtain 



igfi B 2H 7 (e^asAsp 



(1+1) M Q/3 



aP 



ph-R' 



(75) 



SM R = 5M™ + 5Mg +1) = gn B K aP M^_ K , (78) 

with K a p = K^) + K^p X \ The transition rate is based 

upon \Mr + SMr] 2 . Here, |Mr| 2 was taken into account 
above. The interference term can be shown to be pro- 
portional to H z and thus negligibly small. Therefore, 
the field effect is entirely contained in the term |(5A/r| 2 . 
Applying the same procedure as in the previous section, 
one obtains the following addition to the relaxation rate 



aP 



Evaluation of J2 a p \^ap\ 2 for the lowest doublet m 
—S,m' = S yields 



Y,\K a p\ 2 = 2(H 2 z +-H 



ap 



W 2 + A 2 



(80) 



Therefore, the transition rate in the presence of a mag- 
netic field is 



"(H) 



A 2 



W 2 



where T~™* is given by Eq. 16911 and 



6T 



n 3 S 2 Hi ( Wb ) 2 {k B Tf 



no 



6048ft 



(81) 



(82) 



Note that we have used H z -C H± in the last expression. 
One can see from Eq. I|81|l that the relaxation rate due to 
Raman processes dies out when going out of resonance, 
W > A. 



IV. PROCESSES INVOLVING EMISSION OF 
TWO PHONONS 

For processes involving emission of two phonons of, say, 
wave vectors k and q we use the following designations: 

\<t>+) = |«-k,«- q ) , = |"k + l,n q + 1) (83) 

In this case, conservation of energy reads: 



E + — E- = hujh 



ftw c 



(84) 
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The matrix element for this process is, again, the sum 
f the matrix eleme: 
in the second order: 



of the matrix element with Wg? p j, and that with 7ig_ pll 



M E = Mf + M { E 1+1) , (85) 



where according to equations 10 and (1141) . 



x (nk - 1, riq + 1 1^(50^1 n k , n q ) . (86) 

In this case it is convenient to express the phonon matrix 
element as: 



(rik + l,ra q + 1 jc^^l n k ,n q ) 



Q/3 

ph-E 



M, 



/3a 

ph-E 



2M 



Q/3 

ph-E' 



(87) 



where 



M ph-E = ( n q + 1 l^ol n q) ( n k + 1 \8(t>p\ nu) 
h 2 [kxe kXk } a [qxe clXti ] l3 



n q + l)(n k + l). 

(88) 



A. Transitions between eigenstates of S z 

As stated above, the spin-phonon relaxation by emis- 
sion of two phonons may be more important than the re- 
laxation by Raman processes only if the energy difference 
between the spin-states satisfies AE <§; k B T. Provided 
that the energy difference between tunnel-split levels, A, 
is very small, only the relaxation by the emission of two 
phonons between eigenstates of Sz will be considered. 

To this end, we will make use of the spin-Hamiltonian 
(I23[) and the transitions between its eigenstates, \m). 



Adjacent spin levels, m — » m ± 1 



The matrix elements in this case are 



A4 2) = D( T 2m + l)Z mim ±iM*jf_ B 



M 



(i+i) 



= 0, 



with 



'ph-E 

The decay rate is then given by 



ph-E 



(89) 



(90) 



kA k 



T ^ m ±i = ^2 —\M^\ 2 5 (fcj k + hojq + E m±1 - E m )) . 

(91) 



Using the same techniques as in the previous calculations, 
one obtains 



— >m± 



i H 2 m , m ±ilEi[D(2m±l)Y 



ll TT 3 



(k B T) 7 , (92) 



where 



1 f ei , x 3 (e 1 -xfe^ 

I ei = / dx- \, , ; , . 93 

1152 J (e x - l)(e( £ i- x ) - 1 V ' 



2. Non-adjacent spin levels, m — > m ± 2 
In this case, the matrix elements are 



— D ^ m ,m±l^m±l,m±2-Mph-E 



M 



(i+i) 



D' 2 



(±2m + 3)(±2m + 1) i m , m ±i^m±i,m±2 



M. 



±± 



ph-E 



^h-E 



E. 



mil f^m 



Er,-, -f" hiJlc En 



mil J^m 



Em — flLUf 



(94) 



where 

^h-E = \ f 



M xx — , , 

Jw ph-E 1H p h-E 



^-ET'(^E + ^-E) 



The decay rate is 



1 l m.m±l l li±l, m ±2 I E2 D 2 (k B T) 7 



(95) 



, (96) 



with 



1 f" , (e 2 ~x) 3 x 3 e^ 

J-E2 = / dx- 



288 J [e x - l][e^" x ) - 1] 



-(±2m + l)(±2m + 3) 



D/k B T D/k B T 



ei + x ei + e 2 — x 



(97) 



V. DISCUSSION 



Let us now compare transition rates obtained for di- 
rect, or one-phonon, processes H,0J and the rates of the 
same transitions due to Raman processes. 

For spin transitions between adjacent eigenstates of the 
spin Hamiltonian J23J, the ratio of the Raman rate and 
the direct rate is given by 

r „+l {kBT)7[1 _ e -AE/(k B T) ] 

' ■ > I 98 ) 



1 D 



8f{AEf 



where AE = E m ~ E m+1 . In the limit of k B T < AE < 



9 



8d transitions are dominated by direct processes. When 
AE <Si k B T < D one has 



im — >m+l 



D 



\AE 



St 



(99) 



In this case transitions can be easily dominated by the 
Raman processes. For example, for 5=10, m = 5, D — 
0.1K, gf-i B H = IK the energy difference is AE = 2. IK, 
so that at T = 35K the rate of the direct process is 
pOT-vm+i ^ io 5 s _1 , while the rate of the Raman process 
is r^p m+1 ~ 10 6 s _1 . In these estimates we used the 
value of ft ~ 100K. Note that in such molecular magnets 
as Mn-12 and Fe-8 it will be difficult to have the cor- 
responding Raman processes dominant because of large 
distances between adjacent spin levels and small Debye 
temperature. Processes involving the emission of two 
phonons cannot be dominant in any temperature range. 
In some range they can dominate over Raman processes 
but not over direct processes. 

For spin transitions between tunnel-split states of the 



spin Hamiltonian (|44|l at H = 0, the ratio of the Raman 
rate over the direct rate at A <§C k B T is given by 



HO 



960tt 8 



297(25-- l) 2 



k B T 
D 



k B T 
£t 



1 D 

(100) 

Consequently, in zero field and temperatures significantly 
exceeding A, the Raman processes will have much higher 
probability than direct processes. At, e.g., D — 0.5K, 
S = 10, A = 10~ 2 K and T = 5K, the rate of the direct 
process gives T^^ - ~ 10 -4 s -1 , while the rate of the 
Raman process will be ~ 5 • 10 4 s -1 . Note that 

at temperatures where Raman processes dominate over 
direct processes, contribution of the magnetic field to the 
rate, Eq. I)82[l. is small compared to the zero-field rate, 
Eq. PI). 
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